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Abstract. We consider the partition function Z(N;xi, . . . , xp?,yi, . . . ,j/jv) of the square 
ice model with domain wall boundary. We give a simple proof of the symmetry of Z with 
respect to all its variables when the global parameter a of the model is set to the special 
value a = exp(i7r/3). Our proof does not use any determinantal interpretation of Z and can 
be adapted to other situations (for examples to some symmetric ice models). 



1. Introduction 

An alternating sign matrix (ASM) is a square matrix with entries in { — 1, 0, 1} and such that 
in any row and column: the non-zero entries alternate in sign, and their sum is equal to 1. 
Their numbers appear in the so-called Razumov-Stroganov conjecture related to the O(l) loop 
model (0; B 0)- Their enumeration formula was conjectured by Mills, Robbins and Rumsey 
(d), and proved by Zeilberger (fl2l). and almost simultaneously by Kuperberg Qj. Kuperberg 
used a bijection between the ASM's and the states of a statistical square ice model, for which 
he studied and computed the partition function. He also used these tools in (0) to obtain many 
enumeration or cquinumcration results for various classes of symmetries of ASM's, most of 
them having been conjectured by Robbins (fioh . The same method was recently used to obtain 
the enumeration of ASM's invariant (@) or quasi-invariant ([l]) under a quarter-turn rotation. 
A property useful in all these works, which was established by Stroganov (fill), states that 
the partition function Z(N; x%, . . . , xn, yi, ■ ■ ■ , Un) of the (unrestrited) square ice model with 
domain wall boundary is symmetric in all its variables when the global parameter a of the model 
is set to the special value a = exp(i7r/3). We give here another proof of this result. This proof 
is somehow "elementary" since it does not use any determinantal interpretation. Moreover, it 
can be adapted to other cases, for example to some symmetric models, to models in which there 
can be lines that carry more than one spectral parameter, or in which some fixed oriented edges 
are specified. 

This paper is organized as follows: in Section[2j we recall the definitions of the ice models and 
of their associated weights; in Section [31 we give and prove the main result, i.e. the symmetry 
of Z\ in the last section, we present how we can obtain in the same way symmetry properties 
for other ice models. 



2. Definitions 

2.1. Notations. We recall here the main definitions and refer to {!) for details and many 
examples. 
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Let a € C be a global parameter. For any complex number x different from zero, we denote 
x = and we define: 

(f ) a(x) — x — x. 

If G is a tetravalent graph, an ice state of G is an orientation of the edges such that every 
tetravalent vertex has exactly two incoming and two outcoming edges. 

A parameter x ^ is assigned to any tetravalent vertex of the graph G. Then this vertex 
gets a weight, which depends on its orientations, as shown on Figure [TJ 



cr{a 2 ) c(o 2 ) a(ax) cr(ax) a(ax) a{ax) 
1 -1 

Figure 1. The 6 possible orientations, their associated weights, and the cor- 
responding entries in ASM's 

It is sometimes easier to assign parameters, not to each vertex of the graph, but to the lines 
that compose the graph. In this case, the weight of a vertex is defined as: 



V 



xy 



When this convention is used, a parameter explicitly written at a vertex replaces the quotient 
of the parameters of the lines. 

The partition function of a given ice model is then defined as the summation over all its 
states of the product of the weights of the vertices. 

To simplify notations, we will denote by Xjsr the vector of variables (x%, . . . , xjv)- We use 
the notation X\x to denote the vector X without the variable x. 

2.2. ASM's and square ice model. We give in Figure [2] the ice model corresponding to 
(unrestricted) ASM's and its partition functions. The bijection between ASM's and states of the 
square ice model with "domain wall boundary" is now well-known (c/. HI)). The correspondence 
between orientations of the ice model and entries of ASM's is given in Figure [TJ 



x N - 



Z(N;xi, ...,x N ,yi,.. .,y N ) 



x 2 - 

X\- 



2/i V2 yN 
FIGURE 2. Partition function for ASM's of size N 
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2.3. Yang-Baxter equation. To deal with partition functions of ice models, the crucial tool 
is Yang-Baxter equation, that we recall below. 

Lemma 1. [Yang-Baxter equation] If xyz 

x 

(2) 



3. Main result 
The following result has been obtained by Stroganov (fill). 
Theorem 2. (Stroganov) When a = loq = exp(i7r/3) ; the partition function Z(N; Xn,Yn) is 




symmetric in all its variables xi, ■ 



To prove this result Stroganov (jlll ) uses a determinantal interpretation of Z. We want here 
to give a proof that only uses Yang-Baxter equation to study the partition function. 

The method used was introduced by Kuperberg (jH) : observe that Z is a Laurent polynomial, 
then give enough specialization of one of its variable to imply the desired property. 

3.1. Laurent polynomial. Since the weight of any vertex is a Laurent polynomial in the 
variables x^s and y^s, the partition function Z is a Laurent polynomial in these variables. 
Moreover it is a centered Laurent polynomial, i.e. its lowest degree is the opposite of its 
highest degree (called the half- width of the polynomial). Since any row and column of an ASM 
has at least one non-zero entry, which corresponds to a constant cr(a 2 ), we get the following 
property. 

Lemma 3. The partition function Z(N; Xn, Yv) is a Laurent polynomial in any of its variables 
of half-width N — 1 and of parity the parity of N — 1. 

3.2. Partial symmetry. The following lemma gives a (now classical) example of use of the 
Yang-baxter equation. 

Lemma 4. 

(3) 



y- 



Proof. We multiply the left-hand side by a(az), with z = axy. We get 

y- 



o~(az) 



--Z- rx 



X- 

y- 



■3C 

•:xr 
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□ 

As a consequence, we get the "partial symmetry" of Z, which is true whatever the value of 
the global parameter a. 

Lemma 5. The functions Z(N; Xm,Ym) is symmetric separately in the two sets of variables 
Xn and Ypf. 

3.3. Specialization. 

Lemma 6. [specialization of Z; Kuperberg] // we denote 

A(y 1 ,X N \x 1 ,Y N \y 1 ) = J| o-(ax k y~i) J| a(a 2 y 1 y k ), 

2<k<N l<k<N 

A(y 1 ,X N \x 1 ,Y N \y 1 ) = a{ay{x k ) a(a 2 y k y 1 ), 

2<k<N l<k<N 

then we have: 

(4) Z(N;Eyi,X N \ Xl ,Y N ) = A( yi ,X N \xi, Y N \ Vl )Z(N - l;X N \x u Y\ yi ) 

(5) Z{N;a yi) X N \x u Y N ) = A{y x , X N \xi, Y N \ Vl )Z{N - 1; X N \x u Y N \ Vl ). 

Proof. We recall the method to prove equation ([4]). We observe that when X\ = ay±, the 
parameter of the vertex at the crossing of the two lines of parameter x\ and y\ is a. Thus the 
weight of this vertex is a(aa) = a(l) = unless the orientation of this vertex is the second one 
on Figure [1] But this orientation implies the orientation of all vertices in the row x\ and in the 
column yi, as shown on Figure [3] The fixed part gives the partition function Z in size N — 1, 
without parameters X\ and j/i, and the weights of the fixed part gives the factor A{. . .). 



x N - 



■Tl 



x 2 - 
ayi- 



Xi 



ayi- 
x N - 



X2- 



yi 



UN 



y\ 



UN 



Figure 3. Fixed edges for (g]) on the left and © on the right 



The case of ([5]) is similar, after using Lemma [5] to put the line x\ at the top of the grid. 

□ 

3.4. Conclusion. We are now in a position to complete the proof of Theorem [2l From now 
on we set the global parameter to thz special value a — exp(i7r/3). This special value of the 
parameter implies the following equalities: 

(6) a(a) = u(a 2 ) a(a 2 x) = — a{ax) = <j(ax). 
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The proof of Theorem [2] is done by induction on N. We easily check the property for N = 1. 
Now if > 2, because of Lemma we want to obtain the symmetry of Z with respect to 

£2,2/2- 

We use Lemma [3] to reduce the proof of Theorem [2] to the proof of the symmetry of Z in at 
least N (independent) specializations of the variable x±. When x\ = ayi, Lemma [6] gives: 

(7) Z(N; a yi ,X N \ Xll Y N ) = A( yil X N \ Xl ,Y N \ yi )Z(N - 1; X N \ Xl , Y N \ yi ). 

By recurrence, we have that Z{N — I; Xn\xi,Yn\ui) is symmetric in x 2 ,y 2 . The terms of 
A(yi, Xpf\xi, Ypj\yi) involving the variables £2,2/2 are: 

a(ax 2 y\) o-{a 2 yiy 2 ). 

We use equation |6j) to write: 

a(ax 2 yi) cr(a 2 y 1 y 2 ) = v(ax 2 y~i) <r(ay 2 yi) 

which is clearly symmetric in the variables £2,2/2- We obtain in the same way the symmetry 
when x\ = ay\\ and by Lemma[5]we get this symmetry for the 2(A— 1) special values X\ = a^yk 
for k = 1, 3, . . . , N, which is more than enough to conclude the proof of Theorem [2J 

4. Other ice models 

The method used to prove the global symmetry of Z may be adapted to other ice models. 
We illustrate this with half-turn symmetric ASM's (HTASM's). This example shows how our 
method can be used to prove global symmetries, or partial symmetries (for example when a 
line of the ice model carries two different spectral parameters -which breaks the homogeneity 
of the partition function-, or when the orientation of an edge is fixed). 

4.1. HTASM's — notations and results. The ice models corresponding to HTASM's were 
introduced by Kuperberg (|5|) for the even size and by Razumov and Stroganov (|9D for the odd 
size. We recall these models on Figure [4] The spectral parameters on this figure are slightly 
different from the original ones to better suit the proof. The dotted lines mean a change of 
parameter: on one side the parameter is x, whereas on the other side it is y. 

We will prove the following theorem, for the special value a = exp(i7r/3), and for the spe- 
cialization x = y in the even case (which corresponds to the original definition of Zyvt{2N). 

Theorem 7. When the global parameter a is set to the special value — exp(i7r/3), the function 
Zrt(2N; Xn-i, (xn,xn), Yn) is symmetric with respect to all its 2N variables, and the function 
Zht(2A^+1; Xpf, (x, y), Y/v) is symmetric with respect to its 2N variables x%, . . . , xn, 2/i, • ■ ■ , 2/JV • 

The property on Zht(2A^ + 1) may be deduced from the main result of (9;). The assertion 
about Zin(2N) is new. 

4.2. Proofs. Since the method is essentially the same as in the case of unrestricted ASM's, we 
shall give the main steps of the proof and only insist on the difference with the previous case. 

4.2.1. Laurent polynomials. To deal with Laurent polynomials of given parity in the variable y, 
and thus divide by two the number of required specializations of this variable, we shall group 
together the states with a given orientation (indicated as subscripts in the following notations) 
at the edge where the parameters x and y meet. 

So let us consider the partition functions Z^ T (2N; Ajv-i, (x, y)^^) and Z^ T {2N; Xjv-i, {x, y), Yjv), 
respectively parts with the parity of N — 1 and of N of Zm(2N; Xjv— 1, (x,y),Y^) in y; and 
Z^ T (2N+ 1; Xn, Ym, {x, y)) and ^^(2^+1; Xn, Yjv, (x, y)), respectively parts with the parity 
of N and of A - 1 of Z RT (2N + 1; X N , Y N , (x, y)) in y. 
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= Zkt(2N + l;xi, ...,3^,3/1,.. - ,J/JV) (a?> J/)) 



2/1 Vn V 

Figure 4. Partition functions for HTASM's 



Lemma 8. The functions z£ T (2N;X N -i,(x,y),Y N ), Z^ T (2N;X N -i,(x,y),Y N ), Z^ T (2N + 
1; Xn, Yn, {x, y)) and Z^ T (2N + 1; Yjy, (a;, y)) are centered Laurent polynomials in the vari- 
able y, odd or even, of respective half-widths N — 1, N, N — 1, and N . 

4.2.2. Partial symmetries. Since Lemma |4] may be easily adapted to HTASM's, we get: 

Lemma 9. The functions Z^ T (2N; X N -i, (x, y), Y N ), Z^ T (2N;X N -i,(x,y),Y N ), Z^ T (2N + 
1; Xn, Yn, (x, y)) and Z^(2N + 1; Xpf, Yn, (x, y)) are symmetric separately in the variables Xi 
and in the variable yi. 

Moreover if we specialize x = y = xn in the even case, we get the symmetry with respect to 
the set X]\r. 

We now have to deal with the symmetry in the couple (x, y), in the even case. 
The easy transformation 

(8) 5>Q = (a(az) + a(a 2 ))(^ + 

together with Yang-Baxter equation @ gives the following lemma. 

Lemma 10. 

a (a 2 ) + cr(xy) x _ 
V- 



(9) 
(10) 

(11) 



y- 

X- 



v- 

x- 

y- 

X- 



(j{a?yx) 
o-{xy) x- 



a(a 2 yx) y _ 

°jxy) x- 
a(a 2 yx) y _ 



g(a 2 ) x- 

a(a 2 yx) y _ 

o-(a 2 ) x- 

a(a 2 yx) y _ 
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We deduce from this lemma the following property of pseudo-symmetry in (x, y) for the 
functions Z^ T (2N; Xif-i, (x, y),Y N ) and Z^ T (2N; X N _ 1: (x, y),Y N ). 

Lemma 11. For* = 3,5 and □ = 3 respectively, we have 

(12) Z^ T (2N;X N ^,(x,y),Y N ) = \ a(a 2 )Z^ T (2N; X N ^, (y,x),Y N ) 

+ a{xy)Z§ T (2N; X N _ U (y, x), Y N ). 

4.2.3. Specializations. We now give specializations of the functions Zht in the variable x or y. 
Lemma 12. [specialization of Zwr] If we denote 

A 1 H (xi,X N \x 1 ,Y N ) = Yl o-{a 2 xix k ) JJ a(ay k Xi) 

l<k<N l<k<N 

A H (xi,X N \xi,Y N ) = J| a(a 2 x k x 1 ) ]Q er(axiy fe ) 

l<k<N l<k<N 

A° H (yi,X N -i,Y N \yi) = ]J a(ax k y x ) vi^ViVk) 

l<k<N-l l<k<N 

A H {yi,X N ^x,Y N \yx) = ]J[ a(ay{Xk) J| a(a 2 y k y 1 ), 

l<k<N-l l<k<N 

then for * = ~~k *| and □ — 3, 3 respectively, we have 
(13) 

Z^ T (2N + 1; Xn, Yn, (x, axi )) =A 1 H (xx,X N \x u Y N )Z% T (2N; X N \x u (an, a:), Fjv) 
(14) 

(x, axi)) =A ff (an, Xjv\xi, >Ar)ZH T (2A r ; -Xjv\an, (x, xi), Yjy) 

(15) 

Z^ T (2N; X N _ U (x, a yi ), Fjv) =a(axy 1 )A° H {y 1 ,X N _ 1 ,Y N \y 1 )Z° T {2N - 1; Xjv_i, 1^2/1, {x, y x )) 
(16) 

Z° T (27V; X iV _ 1 , (ayi, y), Y N ) =o-(a yi y)A H (y 1: X N ^,Y N \ yi )Z^ T {2N - 1; X N - X ,Y N \ Vl , {y, yi )) 

Proof. The method is almost the same as the one used to prove Lemma [5] An extended proof 
may be found in (jl|). □ 

4.2.4. Conclusion. We are now in a position to conclude the proof of Theorem[7] From now on 
we set the global parameter a to the special value a = exp(i7r/3). 

Lemma fallows us to reduce the proof of the theorem to the proof of the following assertion: 
the functions Z^ T {2N; X N _ lt {x, y), Y N ), Z^ T {2N; X N _ U {x, y),Y N ), Z^ T {2N+l;X N) Y N ,(x,y)) 
and Z^ T (2N + l;Xpj, Y/v, (x, y)) are symmetric with respect to x 2 ,y 2 . 

The proof is done by induction on N. We easily check the property for small values of N. 
We suppose the property true in size 2N, and we consider the odd case 2N + 1. 

Because of Lemma |8l we have to check the property for enough specializations of the variable 
y. We use Lemma [T2l equations (fT3|) and (fl"3)) give us the specializations y = a ±1 xi. As in the 
case of Z, the recurrence gives us the symmetry of the factor Zyvt{2N), and the symmetry of 
the factor A is settled through |6]) which implies: 

a(ax 2 yi) cr(a 2 y 1 y 2 ) = a(ax 2 yi) cr(ay 2 yi) 

the right-hand being clearly symmetric in x 2 ,y 2 . 
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Using Lemma [9] we get 2(N — 1) specializations satisfying the symmetry, which is enough to 
imply the symmetry in full generality as soon as N > 2. 

Now we use the property in size 2N + 1 to prove it in size 2N + 2. The method is the 
same, with the only difference that we have to use Lemma [TU] because equation (jT6j) gives a 
specialization for the variable x. 

The proof of Theorem [7] is now complete. 

Aknowledgment. The author would like to thank P. Duchon for valuable discussions and 
for the use of his figures of ice models. 
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